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Abstract: Resummation of hadron collision cross sections, when the measurement im¬ 
poses a hierarchy of scales, relies on factorization. Cancellation of Glauber / Coulomb 
gluons is a necessary condition for factorization. For Drell-Yan-like processes, the known 
proofs of cancellation of Glauber gluons are not applicable when jet vetoes are introduced, 
via jet algorithms or event shape variables such as the beam thrust. A priori, this does 
not rule out the possibility that an unknown new cancellation mechanism exists, or the 
possibility that a generalized factorization formalism is correct. To resolve the questions, 
we construct a direct counter-example in QCD with scalar quarks, contradicting any form 
of factorization in which the two collinear sectors are decoupled from each other. In the 
counter-example, decoupling of the two collinear sectors implies zero dependence of the 
beam thrust distribution on the longitudinal spin of the incoming hadrons, but we find 
a non-zero spin asymmetry at leading power due to Glauber gluons exchanged between 
spectators. We discuss implications for resumming large logarithms from jet vetoes. 
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1 Introduction 

For the Drell-Yan (DY) process, the cancellation of Glauber gluons was a major difficulty 
in proving factorization [1—3]. In QCD factorization at leading power, by the Collins- 
Soper-Sterman (CSS) formalism [4-6] or soft collinear effective theory (SCET) [7-10], soft 
gluons decouple from the dynamics of collinear particles through the eikonal approximation. 
However, the eikonal approximation is not applicable to soft gluons whose momenta are 
dominated by transverse components, called Glauber / Coulomb gluons. This issue is of 
direct relevance for the resummation of jet veto logarithms at hadron colliders [11-16], 
which requires re-factorization of the Drell-Yan cross section when the jet veto scale is 
much lower than the hard scale. 

In the simpler cases of the Sudakov form factor and semi-inclusive deeply inelastic 
scattering (SIDIS), cancellation of Glauber gluons is achieved at the amplitude level, by 
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deforming integration contours [17] away from the Glauber region. In such cases, if the 
Glauber region is included in the calculation, the cancellation happens between the Glauber 
region and the subtraction of overlap between the Glauber region and other regions [6, 18, 
19]. Such cancellation is by no means automatic; it places strong constraints [6, 20] on 
the choice of Wilson lines in the definition of the soft and collinear functions, including 
the directions (past-pointing / incoming versus future-pointing / outgoing) and possible 
rapidity regularization. 1 

For the Drell-Yan process, the cancellation is more involved, due to the presence of 
both initial-state and final-state poles, in both collinear sectors. Incoming Wilson lines are 
chosen, to be compatible with contour deformation away from initial-state poles (which 
is responsible for the sign flip of the Sivers function between DY and SIDIS [23]), while 
final-state poles that obstruct contour deformation are canceled after summing over cuts 
and integrating over certain momentum components. While there are earlier proofs of 
cancellation of Glauber gluons, the CSS proof in [3], based on both plus- and minus- 
lightcone ordered perturbation theory (LCOPT), is the most powerful one, because its 
applicability extends beyond leading-twist massless parton scattering [24], with important 
phenomenological applications to, e.g. subleading-twist quarkonium production [25]. 

The CSS proof originally required integrating over the transverse momenta of the 
partons, but it was subsequently realized that in the Feynman gauge, the proof carries 
through [6, 26] for transverse momentum dependent (TMD) factorization for the Drell- 
Yan process. However, for factorization of “isolated” Drell-Yan production with measured 
hadronic event shape variables such as transverse energy [13, 27, 28] and beam thrust [29- 
31] (see also [32]), the existing proofs are not applicable, as shown in [33] which explored 
connections with multi-parton interactions. 

In fact, one of the crucial last steps of the CSS proof is integrating over the virtuality of 
the active partons (after summing over cuts in LCOPT). This step is directly broken by a 
measurement of the beam thrust variable, as the factorization proposed in [29] involves the 
virtuality-dependent PDF, also called the beam function. To address the miscancellation 
of Glauber gluons, Ref. [34] introduced new jet veto observables that are designed to be 
less sensitive to such factorization-violating effects. 

We will borrow the terminology “generalized factorization” for hadron-hadron colli¬ 
sions proposed in a slightly different context, TMD factorization. Generalized factorization 
is to be distinguished from “standard factorization” for hadron-hadron collisions, the lat¬ 
ter of which assumes cancellation of Glauber gluons and always defines soft and collinear 
functions using past-pointing / incoming Wilson lines carrying the color charges of the 
active partons. For example, the factorization of beam thrust in [29] using SCET should 
be characterized as standard factorization." Generalized factorization, in a narrow sense, 

lr The compatibility between contour deformation and rapidity regularization by off-lightcone Wilson lines 
[6] has been studied in the aforementioned references, but it should be possible to extend the studies to 
other rapidity regulators [21, 22] used in the SCET literature, which bear more resemblance to dimensional 
regularization. 

2 The equivalence of leading-power soft collinear factorization derived from traditional QCD methods 
and from SCET, for sufficiently inclusive observables that guarantee the cancellation of Glauber gluons, 
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involves modification of Wilson lines in collinear and soft functions [40-43], but in a gen¬ 
eral sense, can be any factorization with a sensible spin structure [44]/’ In the model field 
theory considered in this paper, the colored active partons and the produced heavy uncol¬ 
ored particles are scalars carrying no spin indices, and generalized factorization is liberally 
defined as any factorization of the cross section into a product / convolution of scalar hard, 
collinear, anti-collinear, and possibly soft functions. We will make the definition precise in 
Section 3.1. 

To disprove a leading-power factorization statement, it is sufficient to demonstrate, at 
some fixed order in a s , that the leading-power part of the cross section contradicts the 
prediction from factorization. Since QCD factorization relies on arguments that are appli¬ 
cable to all unbroken gauge theories, irrespective of e.g. gauge groups and matter contents, 
it is sufficient to find a contradiction in a model field theory that allows easy calculation. 
Model field theories involving polarized scattering, again in the slightly different context 
of TMD factorization, have been used to show the violation of both standard factorization 
and generalized factorization for hadron production at small transverse momentum [44-47]. 
In particular, Ref. [44] exploited discrete symmetries to study the violation of factorization, 
and this approach will be adopted in our study. 

In this paper, we will consider Drell-Yan-like scattering in a model field theory, and 
study spin asymmetries in the doubly-differential beam thrust distribution. For this special 
model and observable, the vast majority of diagrams vanish, allowing a clean calculation 
of factorization-violating effects. The goal of the study is two fold. First, we would like 
to give an explicit demonstration that standard factorization is violated, which was shown 
by [33] to be extremely likely. Second, we would like to show that it is not possible for 
a generalized factorization theorem to hold, since decoupling of the two collinear sectors 
necessarily leads to zero spin asymmetry, while the calculations in this paper find a non¬ 
zero spin asymmetry. The non-cancellation of Glauber gluons found in this paper only 
happens above the jet veto scale, so collinear factorization is still valid if the jet veto scale 
is perturbative, but our ability to resum large logarithms in the hard scattering function 
will be compromised. 

Ref. [34] introduced new jet algorithm-based observables, such as the “jet beam thrust”, 
which are designed to be less sensitive to Glauber effects, while preserving the rapidity- 
dependent nature of the beam thrust variable. Other recent research [44-52] investigated 
the violation of QCD factorization in contexts other than, or wider than, the Drell-Yan pro¬ 
cess. For perturbative resummation, the logarithmic order at which factorization-violating 
effects start has been discussed for top quark pair production at low transverse momentum 
[53, 54] and dijet event shapes [55]. 

This paper is organized as follows. In Section 2 we give a description of the model field 
theory and the observable measured in our thought experiment. In Section 3 we explain why 
a non-zero spin asymmetry would contradict both standard factorization and generalized 
factorization. In Section 4 we check that up to O(a^), factorizable diagrams, as well as the 

has been demonstrated extensively, for example in [35-39] 

‘Tor example, the large component of the Dirac spinor of a collinear quark should be projected out, and 
gluons that enter the hard scattering should carry two possible transverse polarizations. 
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Figure 1 . The leading-order Feynman diagram for 7 + 7 —► $ + A' in the model theory. Here 
X = ()> + (jf. 


vast majority non-factorizable diagrams, do not contribute, due to cancellations associated 
with this special observable. In Section 5 we evaluate the only 0(a 2 s ) cut diagram left, 
which involves one Glauber gluon exchanged on each side of the cut. The evaluation gives a 
non-zero spin asymmetry, which is the desired result. Some discussions are given in Section 

6 . 

2 The model for showing factorization breaking 

2.1 Model field theory 

We consider QCD with a hypothetical massless complex scalar quark cf (instead of a Dirac 
fermion quark as in real QCD) under the fundamental representation of SU(3). The scalar 
quark also carries an electric charge of +e. The produced heavy particle (analogous to 
7 * and Z° in the real Drell-Yan process) is a neutral color-singlet scalar with mass M, and 
couples to the light scalar quarks via the interaction Lagrangian —g§&cfi* <f>. The photon 7 
serves the role as a ’’hadron” which can split into an active scalar quark that participates 
in hard scattering, and a spectator scalar quark going down the beam pipe. We consider 
the Drell-Yan process 7 + 7 —> + X. The lowest-order diagram is 7 + 7 —> <f> + 0 + <)>*, 
shown in Fig. 1. 

2.2 Observable - spin asymmetry in doubly differential beam thrust 

We will use the beam thrust variable [29, 56] 

TB = J^Yl\ pT i\ e ~ lVi ~ vl , ( 2 - 1 ) 

i 

where y and M are the rapidity and invariant mass of the Drell-Yan pair (actually a heavy 
scalar in our model), and the index i runs over every detected hadron. Requiring rg < 1 
strongly restricts hadronic activity, especially in the central rapidity region. It is useful to 
consider the doubly-differential distribution [29] in (tr,tx), where tr receives contribution 
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from only the right hemisphere (y* > y) and t r receives contribution from only the left 
hemisphere (y* < y). 

Even though the arguments in this paper also apply to jet algorithm-based vetoes, 
for concreteness we will consider the doubly-differential beam thrust distribution from 
the scattering of two photons. Our thought experiment involves incoming photons with 
several polarization configurations, from which we obtain the double spin asymmetry in the 
distribution. We will do the calculation for a heavy scalar $ produced at rapidity y = 0, 
and with the left and right hemispheres having the same beam thrust r R = tl = 2 t r <C 1. 
The phase space integral for the square of the amplitude in Fig. 1 is, in the approximation 
that k\ is plus-collinear and k 2 is minus-collinear, 
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dT R dr L dy 
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The delta functions in Eq. (2.2) force 
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(2.3) 

(2.4) 


d (\kir \ 2 - |&ot| 2 ) d [\k 2 T \ 2 — |& ot | 2 ) \M\ 2 . (2.5) 


For the scattering of two photons, we define the absolute double spin asymmetry as 


°asym — (off + <T 4t ~ a ft ~ a li) /(2-6) 

where up and down arrows denote right and left polarizations. The relative double spin 
asymmetry is defined as the above expression divided by the unpolarized cross section 


°unpoi — (044 + °"4t + “I" Tld) /^- (2-7) 

When a photon traveling in the 2 direction is right polarized the polarization vector is = 
(e°,e*,e|,e|) = (0,1, i, 0)/\/2. So efe* u = (-g^r _ * e nr«*r) / 2 , where e xy = -e yx = 1. 
Similarly, for a left-polarized photon, = ( y —g yTl " T + ie yTUT ) /2. The half difference 

between the polarization sums is 

\ (eft" - 44“) = ( 2 . 8 ) 
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To obtain the (absolute) double spin asymmetry in the distribution, we can replace \M\ 2 
in Eq. (2.5) by 

IMlLym = \e^ TVlT E^ TV2T M^^M^^, (2.9) 

while to obtain the unpolarized cross section, we use the averaged squared matrix element 

l-^lunpol = \g^ T u 1T 9^ T u 2T M^^M^^, (2.10) 


3 Outline of the proof by contradiction 


3.1 Definition of standard factorization 


By standard factorization, we mean factorization derived by assuming the cancellation 
of spectator-spectator Glauber gluon exchanges. For the doubly differential beam thrust 
distribution (in the left and right hemispheres) in Drell-Yan production of a system of 
invariant mass M, at hadronic c.o.m. collision energy E, the factorization formula in double 
Laplace moment space is [29, 31] 


roo roc 

/ dr R e~ NR TR / cLtl e~ Nlj TL 

Jo Jo 

{ M 2 


d 3 a (H 1 + H 2 ->$ + X) 


dT R dr L dy 


y =o 


. M 2 , \ ~ ( M 2 . \ ~ ( M 2 

= H y-jj-.a, ( H )j Bl j B 2 

- ( M 2 M 2 .A ^ / 1 1 \ 

x {w? J ’W^' a,M ) + \Ni’Ifl)’ 


, X, Ot S (/-i/ ) 


(3.1) 


where x = M/E is the Bjorken variable for both collinear sectors, as we imposed y = 0. 
H is the hard function. B\ and B 2 are the two moment-space collinear functions, also 
called beam functions in the literature, for the collinear sectors initiated by the incoming 
hadrons H\ and H 2 , respectively. S is the moment-space soft function. The beam function 
is the “virtuality-dependent PDF” [29], originally defined using SCET fields with implicit 
zero-bin subtraction [57]. In Appendix A we define the beam function using the scalar 
QCD fields used in this paper. 


3.2 Definition of generalized factorization 

We first re-write Eq. (3.1), dropping explicit dependence on y,f and a s (yf), and combining 
H and S into a function R, 


/ dr R e NrTr / dr^e 

/ o Jo 


- Nl Tl dJa (H 1 + H 2 ^<^ + X ) 
dT R dr L dy 


y =o 


= H ( M 2 ) (M 2 , N r , x) B 2 (M 2 , N l , x) S (M 2 , N r , N l ) + O 


1 1 

Nr’Nl 


Bl (M 2 , N r , x ) B 2 ( M 2 ,N l , x) R (M 2 , N R , N L )+o(^-,^-Y (3.2) 
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Generalized factorization is defined by Eq. (3.2), with the following requirements. B\, 
B 2 , and R are allowed to be arbitrary functions that may differ from how they are defined 
in standard factorization. 1 The only dependence on the incoming state Hi, including 
the species of the hadron and the polarization, should be contained in B i, and the same 
condition is imposed on H 2 and i? 2 - In other words, we require the two collinear sectors to 
be decoupled. The various functions on the R.H.S. of Eq. (3.2) generally still depend on 
additional auxiliary variables that are not shown, including the factorization scale Bf and 
the directions of Wilson lines involved in the definition of these functions. 

It is the aim of the next subsection to show that generalized factorization is incom¬ 
patible with the result of this paper; standard factorization, a special case of generalized 
factorization, is thus also violated. Before doing so, we give an example of generalized fac¬ 
torization that differs from standard factorization. In the standard factorization formula, 
Eq. (3.1), the soft function is defined using incoming Wilson lines [29]. As shown in [32], 
up to O(a^), if outgoing Wilson lines are used instead, the soft function remains the same 
up to By flipping the directions of Wilson lines in the soft function, we turn Eq. 

(3.1) into a generalized factorization formula which coincides with standard factorization 
at the first few a s orders, but differs at higher a s orders. 

3.3 Violation of generalized factorization 

Under the generalized factorization formula Eq. (3.2), the double longitudinal spin asym¬ 
metry in the factorized beam thrust distribution is, in schematic moment-space factorized 
form, 

b\) ( B\ - U,j) R, (3.3) 

while the corresponding expression for single longitudinal spin asymmetry is 

^h b\ - i?2 npo1 ^- (3.4) 

In Eqs. (3.3) and (3.4), we used “f”, and “unpol” to denote right polarization, left 
polarization, and no polarization, respectively. It is an immediate consequence of parity 
conservation of our model theory that the single spin asymmetry in Eq. (3.4) must be 
vanishing. Since the unpolarized scattering cross section is non-zero, the unpolarized beam 
function B l, 111 * 301 in Eq. (3.4) cannot be vanishing (except at isolated points, assuming the 
function is analytic). So B[ — must be almost everywhere zero, which means everywhere 
zero if the function is analytic. 0 This implies that the double spin asymmetry, given in 
the schematic factorized form Eq. (3.3), must also vanish at leading power to all orders 
in ot s . But we need to consider non-factorizable contributions from the Glauber region in 
spectator-spectator interaction, and its overlap with other regions (to be subtracted); if 

4 For example, B 1 in Eq. (3.2) can depend on M 2 and Nr independently, while B 1 in Eq. (3.1) is written 
in a form that can only depend on these two variables through the combination M 2 /Nr. 

'Indeed, for the special case of standard factorization, we can check from the definition of the beam 
function that B 1 and B 2 do not depend on the polarizations of the hadrons. 


- 7 - 



the sum of these contributions is non-zero, as we will find in the subsequent sections, we 
obtain a contradiction to generalized factorization. 6 


4 Vanishing diagrams 


4.1 Vanishing LO diagram 

We know that factorizable diagrams give a vanishing contribution to the double spin asym¬ 
metry in our model, but we will explicitly verify that the LO contribution from squaring 
the diagram in Fig. 1 vanishes, to introduce notations and demonstrate cancellations that 
are also applicable to some non-factorizable diagrams. For brevity, we write 

D (p) = -w—■, 
p z + le 

e(p 1 ,p 2 ) = Wm 2 tPi 1 T P2 2T - 

The diagram evaluates to, given that p\ and p 2 have no transverse components, 

= (-i)g^2iek^ T (-2ie)k^ 2T D(k 1 )D(k 2 ) 

= -4 ie 2 gs> fcf' r k^ 2T D (A*) D (k 2 ). 


(4.1) 

(4.2) 

(4.3) 


In this expression we omitted the color factor 5 a b, with a and b being the color indices for 
the scalar / anti-scalar pair in the final state. Using Eq. (2.9), the resulting double spin 
asymmetry in the squared matrix element is, 


^ e /UT V 1 tU*2T V2T-M LQ 


IT M2T 


(M^^r^e(kj,kf)e(klkj)=0. 


(4.4) 


In contrast, the unpolarized spin-summed squared matrix element is, using Eqs. (2.3) for 
zero rapidity, 

\4 


1 


l-^Lolunpol — 2Ca ' a 9/j.it<'it9h2TV2T-Mi j o ^ 2T {-M. 


vitv 2 t)* _ 2C A ■ 4e 4 g| (1 X) 


\kir\ 2 | k 2 T | 2 


(4.5) 

In this expression we include the color factor C A from the final-state color sum, and an 
overall factor of 2 to account for the possibility of reversing the complex scalar arrow (i.e. 
swapping scalar and anti-scalar). The explicit momentum components in Fig. 1 have been 
used to evaluate D(k \) and D(k 2 ). Using Eq. (4.5) as the squared matrix element in the 
phase space integral Eq. (2.5), we obtain the LO unpolarized beam thrust distribution 


d 3 ULO 


M 2 


d,T R dT L dy 
where we defined 

d?k\T f d 2 k 2 T 


y= 0 ,TR=TL=T B /2 


16ttE 2 
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(4.6) 
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\kir\ \k 2 rY 


(47r) 2 \k 0T [ 


(4.7) 


6 We will actually show that the distribution is non-zero at some ( tb, tl), but this is sufficient to imply 
that the distribution cannot be an identically vanishing function in moment space. 
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Figure 2. The one-Glauber exchange diagram in the model field theory. 

4.2 Vanishing one-loop cut diagram 

Consider the diagram Fig. 2 in which the two spectator lines are connected by one gluon 
that is either a Glauber gluon or a normal soft gluon, in interference with the complex 
conjugate of the LO diagram. We will show that the resulting contribution to the spin asym¬ 
metry has vanishing real and imaginary parts. For more general models and observables, 
the contribution is purely imaginary when the gluon has Glauber-like momentum, and 
cancels with the complex conjugate cut diagram [1, 33].' But in our special model theory, 
even the imaginary contribution vanishes. Fig. 2 evaluates to, again noticing pj = pj = 0 
and omitting the color factor Cpdab, 

iM = - ig* j 4 e 2 g 2 s (Af T - l^ T ) (k? T + l^) (2 Pl - 2k 1 + l ) • (2 p 2 - 2k 2 - l) 

D ( l) D (pi -ki + l)D (p 2 -k 2 -l)D (fci -l)D (k 2 + l) (4.8) 

~ -8 ie 2 g<s> g 2 s (pf - k +) (pf — kf) J d 4 l (kf 1T - l^ 1T ) (ktf 2T + l^ 2T ) 

D ( l) D (p! - h + l) D (pa — k 2 — l) D (h — l) D (k 2 + l ), (4.9) 

where we used approximation l + ,l~ <C Q, applicable to both the Glauber and the normal 

soft region, after the sign. Therefore, the interference between the diagram in Fig. 2 
and the complex conjugate of the LO diagram in Fig. 1 is 

c c V2T ( \jV1T V2T\* 

t niT I / lT t M2T V2T JVl l \- /Vt LO ) 

oc J d A le (fcf, l T ) e (fcj, l T ) D ( l ) D (pi - h + l)D {p 2 - k 2 - l) D (A* — l) D (k 2 + l ), 

(4.10) 

where only /y-dependent terms are shown after the proportional sign “oc”. Recall that we 
would like to measure the doubly-differential beam thrust distribution at some ry = tr, 

7 A rare exception is single transverse spin asymmetry, for which factorization can be violated by the 
exchange of only one Glauber gluon, because the resulting imaginary contribution is multiplied by another 
imaginary factor from the Dirac trace with , to give a real contribution [45] . 



Pi - h 
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with the heavy particle <b produced at zero rapidity. We still need to integrate the squared 
matrix element Eq. (4.10) over the phase space of k\ and k 2 with appropriate measurement 
functions. Consider a particular point in the /-integration volume in Eq. (4.10), for example, 
a point with l x 7 ^ 0, l y = 0 without loss of generality. Then we can flip the sign of the 
y component of k\ without changing any terms in the integrand in Eq. (4.10) except for 
flipping the sign of e (/c^,/ r ). Since jet veto observables are azimuthally symmetric and 
do not generate a preferred y-direction for k±, Eq. (4.10) gives a vanishing contribution to 
the doubly differential beam thrust distribution. At some general value of /, the needed 
change of variable is 

ix uv __ yyux 

kj -+ r it o kf = (kf,k\) - 2 n , ( 4 . 11 ) 

denoting a reflection of k\ in the line through the origin in the ±It direction. The squared 
matrix element in Eq. (4.10) has odd parity under this transformation, while the beam 
thrust variables (tr,tl) are invariant. Therefore the contribution vanishes after phase 
space integration. 

4.3 Vanishing two-loop cut diagrams 

Consider any cut diagram whose lower spectator line has only one gluon attachment, with 
the gluon being either soft or Glauber-like, such as the diagrams shown in Fig. 3. As is 
the case for the one-gluon diagram in Fig. 2, the lower collinear sector only depends on 
the following three momenta, pi, k±, and /. It is also easily checked that, again, at leading 
power, the only numerator factors that depends on It (or the other unlabeled transverse 
loop momenta) are the photon-scalar vertices. 

Therefore, exactly the same transformation as in Sec. 4.1, Eq. (4.11), reverses the sign 
of the cut amplitude, and proves that the contribution is zero after phase space integration 
for k\ and k 2 - 

To prove that the cut diagram in Fig. 4, i.e. a 2-loop diagram with a box in interference 
with the complex conjugate of the LO diagram, vanishes, we need a little more work. Recall 
that a “pinch” in the Glauber region arises when both the active parton and spectator 
lines depend on the Glauber-like exchanged momentum. Since the active quark line in 
Fig. 4 depends on /, but not / 1 , only the overall exchanged momentum l can be pinched in 
the Glauber region ~ (A 2 , A 2 , A), to produce a potentially non-factorizable contribution. 
Meanwhile, it can be checked by IR power counting that the individual momenta l\ and l—l\ 
can be both Glauber-like, or both soft ~ (A, A, A), for the diagram to give a leading-power 
contribution. 8 

As is the case for diagram Fig. 2, the only leading-power dependence of numerator 
factors on It is from the photon-scalar vertices. Here these vertices only depend on It but 

8 If l is Glauber-like, h and l — l\ can also be both plus-collinear or both minus-collinear for the diagram to 
contribute at leading power. But it is not necessary to consider this situation, because the sum of diagrams 
involving a secondary hard vertex is suppressed by Ward identities, as shown by Labastida and Sterman 

[ 58 ]. 
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Figure 3. Example cut diagrams, each with only one gluon attached to the spectator line on the 
lower half of the graph, while more than one gluons may attach to the upper spectator line. 
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Figure 4. Two Glauber gluons exchanged on the same side of the cut, in the model field theory. 
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not 1\t • So the leading-power cut-amplitude is, omitting the color factor Cp and other 
constant factors, 


c c X/t^lT H2T ( \ytVlT V2T\* 

fc MlT I'lT t M2T V2T JVl 2 V^LO ) 

oc J dHe (kf, l T ) e (fcf, l T ) D (A* - l) D (, k 2 + l)D (pi - k x + l) D (p 2 - k 2 - l) 
f d\ D (h) D (l — h) D (p! -h+ h) D (pa - k 2 - h) , 

(4.12) 


where we have collected all the Zi-dependent terms in the third line. 

We briefly comment on the “Glauber-II” region discussed in [29], which has mo¬ 
mentum scaling (A 4 , A 4 , A 2 ) instead of the usual Glauber scaling (A 2 , A 2 , A). The factor 
e e (k%,l T ) in Eq. (4.12) and Eq. (4.10), absent in unpolarized scattering, gives a 

suppression when It is smaller than the usual Glauber transverse momentum. Therefore, 
even if the “Glauber-II” region is relevant for leading-power unpolarized scattering, it is 
not relevant here. Also there will be no singularity from It —> 0, which is important be¬ 
cause the next step is analyzing the cut diagram at fixed It, assuming that the subsequent 
integration over It causes no complication. 

With l being Glauber-like and l\ being either Glauber-like or soft, the only propagators 
that have leading-power dependence on l~ are the two lines immediately connected to the 
lower incoming hadron pi, and the only propagators that have leading-power dependence 
on l + are the two lines immediately connected to the upper incoming hadron p 2 . So at any 
fixed It, we can perform the l + and l~ integrals by contour integration [1, 6], picking up 
the poles by cutting the lines p 2 — k 2 — l and p\ — k\ +1 in Fig. 4, producing an imaginary 
contribution multiplied by a one-loop box (sub-) diagram initiated by the on-shell lines 
p 2 — k 2 — l and p\ — k\ + l. To obtain a real contribution that is not canceled by the complex 
conjugate cut diagram, we need another imaginary contribution from applying cutkosky’s 
rules to the one-loop box diagram, with l\ in the Glauber region. 9 The approximations are 


D{h) 

D(l — li) 


D(pi -h + h) 


D(p 2 -k 2 - h) 


i 

— |^it | 2 + *o 

i 

— | It — Iit\ 2 + *0 

i 

(1 - x)Q ( ^l\ Q + ^l) ~ l^i t ~ hr \ 2 + *0 
i 

(1 - x)Q ((yzyj Q -ll') ~ 1^2 t - 1\t \ 2 + io 


(4.13) 

(4.14) 

(4.15) 

(4.16) 


Following Cutkosky’s rules, we replace D{p\ — k\ + l\) and D{p 2 — k 2 — l\) by corre¬ 
sponding delta functions, and integrate over and if, producing the constant factor 


9 As shown in [19], applying cutkosky’s rules to the direct box diagram alone yields the same imaginary 
part as applying the method of regions to the sum of the direct box and the cross box diagram. 
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Figure 5. The cut diagram with one Glauber gluon exchanged on either side of the cut, i.e. the 
square of Fig. 2. 

(27tz) 2 / [(1 — x)Q] 2 . Now the last line of Eq. (4.12) becomes just the two-dimensional 
f d 2 l\T integral over D(l\) and D(l — l\). At this point, we can re-use the strategy for 
showing the vanishing of the interference between Fig. 2 and the LO graph, and replace 
kj in the first line of Eq. (4.12) by its mirror image in the line through the origin in the 
direction of ±l T , as in Eq. (4.11). The integrand reverses sign, and since this transforma¬ 
tion does not affect the doubly differential beam thrust distribution, the final contribution 
from the cut diagram Fig. 4 is zero. 

5 Non-zero two-Glauber diagram 

We have excluded all cut diagrams whose lower (or upper) spectator line is attached by 
only one soft gluon (normal soft gluon or Glauber gluon) in Section 4.3. We have also 
excluded the diagram involving two Glauber gluons on the same side of the cut, including 
the diagram with a box, Fig. 4, and an unshown diagram with a cross box. This leaves us, 
at the order a 2 relative to LO, the only possible contribution from the cut diagram in Fig. 
5, showing the interference of Fig. 2 with its own complex conjugate. Out of the two soft 
gluons exchanged on the two sides of the cut, one has to be a Glauber gluon to produce 
factorization-violating effects, and then the other one also has to be a Glauber gluon to 
produce a real contribution that is not canceled by the complex conjugate cut diagram 
[1, 33], 

5.1 Reducing to 2D integrals by contour integration 

We show that the cut diagram Fig. 5, from squaring the amplitude in Fig. 2, gives a 
non-zero contribution to the double longitudinal spin asymmetry. 
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As above, we fix the heavy particle $ to have zero rapidity, so that in Fig. 5, x\ = 
X 2 = x, kf = k^ = xQ, and the amplitude Eq. (4.9) can be re-written as 


2 2 /-i \2/-i2 / dl + dl d lx 


= —A.ie 2 g$g 2 s (\ — x) 2 Q 


(2tt) 


(fcf T - P 1T ) (jfc£ 2T + Z M2T ) 


D (0 D (pi - Ail + 0 £> (p 2 -k 2 -l)D {ki — l) D (k 2 + l ), 


(5.1) 


In this expression, with l lying in the Glauber region, the relevant leading-power approxi¬ 
mations are 

i 


D(l) 

D(pi -h+l) 
D(k i - l) 
D(p 2 — k 2 — l) 
D{k 2 - l ) 


—tL + i0' 


(1- 

x)Q | 

f IUt| 2 

{(l-x)Q+ l J 
i 

1 

IF 

4 

1 

£T 

(O 

+ 

css. 

o 

xQ \ 

f lUri 2 

{ (l -X)Q ) 

1 — \k\T ~ It \ 2 + *0 
i 

(1- 

x)Q | 

f |^2t| 2 

i 

;' + ) 

1 — \k 2 T — h\ 2 + *0 


xQ ( ~ (t k ll\ Q - l + ) - \h T ~ h\ 2 + i0 


(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 


We first integrate over l + and l using contour integration, picking up the poles from the 
vanishing of D(p\ — k\ + l) and D(p 2 — k 2 — l). We are left with the It integral, 

' d 2 l T 


iMi 1Tll2T = 4:e 2 g$ g 2 (1 — x) 2 

1 1 


(2vr) 2 

1 


(k£ 1T -V tlT ) (k% 2T + l^ T ) 


\h\ 2 \kiT — h\ 2 \k 2 T + h\ 2 


(5.7) 


So the asymmetry from the cut diagram Fig. 5, using the spin sum formula Eq. (2.9), is 

j2 ' ■ ' d 2 l B 1 


= 2C F C A f^f 


(9^-)2 *'lT c A t 2T V 2T 

27 . r 


j^HiT X2T M2T Y 


T 2 r< „4/i „\4 f d Ia f d Ib 


= 2CpC A ■ 4e 4 g$ g^(l — xy 
where, using the notation in Eq. (4.2), 


( 2 vr ) 2 J ( 2 vr ) 2 asym 


Zasym{klT,k 2 T,l A ,l B ) , (5.8) 


T^asym (&1 Ti &2T> lA, ^b) — £ (&1 — — ^b) £ (k 2 + lA, &2 + ^b) 


Ka| 2 | k±T — l A \ 2 \k 2 T + Ia\ 2 \Ib\ 2 \klT — Zs| 2 \k-2T + Ib\ 2 


(5.9) 


In Eq. (5.8), the previously ignored color factor C\C A is shown, and an overall factor of 2 
is present to account for the possibility of reversing the complex scalar arrow (i.e. swapping 
scalar and anti-scalar) in Fig. 2. Using Eq. (5.8) as |A4| 2 in Eq. (2.5), we obtain 


d 3 




asym 


dr R dT L dy 


M 2 


J/=0,T fl =r i =r fl /2 


167T.E 2 


2C 2 F C A -^g^ 2 gt{l-xfl i 


asym 5 


(5.10) 
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where we defined I asym as the loop and phase space integral over the integrand, Eq. (5.9), 

r _ f d 2 l A f d 2 l B f d 2 kir f d 2 k 2T 

asym J (2tt) 2 J (2tt) 2 J (2tt) 2 J (2tt) 2 

d (\kir\ 2 — \koT\ 2 ^) 6 ^|£;2t| 2 ~ |&0t| 2 ) ^asym {klT, k 2 T, Ia, Ib) ■ (5-11) 

5.2 Cancellation of IR divergences and numerical evaluation 

Consider the l a and l B integrals in Eq. (5.11). Possible IR divergences may arise from the 
vanishing of any of the denominators in Eq. (5.9), if the singularities are not integrable. 
The singularities at l A)B = k\r and 1 Ai b = — k 2 x are integrable, because the numerator 
factors on the 2nd line of Eq. (5.9) vanish at these points. However, there seems to be 
logarithmic singularities when l A or Ib approaches zero. We show that this is not the case, 
because there is a linear suppression when l A or Ib , becomes small, turning the singularities 
into integrable ones. The subsequent argument will be similar to Section 4.1, where we 
used the reflection of k\ to show that the one-Glauber diagram vanishes. This suppression 
also implies that there is no leading power contribution from any regions with much smaller 
transverse momenta, such as the ultra-soft region and the “Glauber-II” region discussed in 
[29], or from the overlap between these regions and the regular Glauber region. 

We use R p o q, defined in Eq. (4.11), to denote the reflection of the two-vector q with 
respect to to the line through the origin in the direction of the two-vector p. We then 
re-write Eq. (5.11) as 


Ttsym — 0 


d 2 kir f d 2 k 2 T 


(2vr) 2 

d?l A 

(2vr) 2 


( 2 7 t ) 2 

d 2 ls 
(2vr 


<5 (\kiT\ 2 ~ |&ot| 2 ) d [\k 2 T\ 2 — I^otT 

2 [-^asym (■ klT,k 2 T,lA,lB ) + -^asym ( Ri a ° kiT, k 2 T, Ia, Ib) ] ■ 


(5.12) 


In this form, the last line can be readily checked to vanish when Ib= 0, and linearly sup¬ 
pressed when Ib is small. The Jacobian factor from the reflection is 1, so we simply need 
to put a factor of 1/2 at the start of Eq. (5.12). We go one step further by reflecting k 2 T 
with respect to the line through l B , recasting Eq. (5.11) into 



Tasym (&1T, Rl B ° &2T, l A, h) + T asym ( Ri a o k\T, Rl B ° k 2 T, lA, Ib ) ] , (5.13) 


where the sum inside the square bracket receive a linear suppression when either l A or l B 
become small, and a quadratic suppression when both l A and l B are made small simultane¬ 
ously. This makes both the points l A , l B = 0 integrable singularities despite the quadratic 
denominators in Eq. (5.9). So the expression Eq. (5.13) is IR finite. UV finiteness is also 
clear by power counting. With both IR and UV divergences absent, Eq. (5.13) can be eval¬ 
uated by straightforward Monte Carlo integration without regularization or subtraction. 
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Using the Vegas algorithm implemented by the CUBA library [59], with 4.2 million points 
sampled, we obtain 

/ asym = (1.58 ± 0.02) ] (5.14) 

( 47 t ) 4 | fe 0T | 

Dividing the absolute asymmetry given by Eqs. (5.10) and (5.14) by the LO unpolarized 
differential cross section given by Eqs. (4.6) and (4.7), we obtain the relative spin asym¬ 
metry 


/ d^asym , d 3 a LO 
\dT R dT L dy dT R dr L dy 


y=0,T R =r L =T B /2 


cWs Gsyrn / Ax) 

(1.58 ±0.02 )C 2 F a 2 s . 


(5.15) 


We have shown that the asymmetry is non-zero, proving that Glauber gluons break factor¬ 
ization for the double spin asymmetry in the doubly differential beam thrust distribution. 


6 Discussion 

We have performed a calculation of factorization-violating effects in the beam thrust dis¬ 
tribution from Drell-Yan-like scattering in a simple model field theory. Any factorization 
in the limit of small beam thrust (corresponding to a stringent jet veto), standard or gen¬ 
eralized, would predict a vanishing double longitudinal spin asymmetry, due to the scalar 
nature of the active quarks in this parity-conserving model. The non-zero result found in 
our calculation is in contradiction to any generalized factorization that separates beam- 
thrust dependence into universal functions. The non-factoring contribution, Eq. (5.15), is 
in fact infrared safe, which shows that collinear factorization is respected to this order, but 
with factorization scale y R = O(^tbs). Logarithms of beam thrust are thus contained in 
the hard-scattering function of collinear factorization, and standard resummation methods 
do not reply. 

By looking at double spin asymmetry in a theory with scalar quarks, the calculation is 
simplified enormously. For example, the diagrams in Fig. 3 and 4 all vanish, and the only 
non-zero diagram Fig. 5 is IR-finite in the Glauber region without regularization. Neither 
of these simplifications hold for unpolarized scattering in real QCD. Since the factorization 
theorem in [29] for the beam thrust distribution is proposed in a manner that applies to 
both unpolarized and polarized scattering in any unbroken gauge theory, our result is a 
counter-example which shows the aforementioned factorization theorem cannot hold true 
in every situation. But strictly speaking, we do not exclude the small possiblity that the 
factorization theorem survives in unpolarized scattering; a dedicated, more complicated 
calculation for unpolarized scattering would be needed to conclusively settle the question. 

The breakdown of generalized factorization would eventually lead to corrections to the 
existing predictions of jet veto resummation calculations [11—16]. The question remains, 
“At which logarithmic order do such corrections start?’' The lowest-order factorization- 
violating diagram in this study involves two spectator lines and two virtual Glauber gluons, 
producing a non-zero result. This would be of order af if we were studying massless parton 
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scattering instead of photon-photon scattering. This result contains no large logarithms for 
two reasons. First, the intrinsic virtuality of collinear particles, anti-collinear particles, and 
Glauber gluons are all of the order Q 2 tb [29]. Second, the order of the diagram is too low 
to acquire Regge-type rapidity logarithms, which will show up at higher orders in ladder- 
type diagrams. But we still need to multiply the result by the hard function (which always 
factorizes, though collinear and soft functions can be entangled by Glauber gluons), with 
double logarithms ~ a™ ln 2n tb due to running from the scale Q 2 to Q 2 tb■ We end up with 
a" +4 ln 2n TB. In unpolarized scattering, there can also be a non-zero contribution when 
two soft / Glauber gluons are exchanged on the same side of the cut, which potentially 
gives one power of a Regge-type logarithm, resulting in a” +4 ln 2n+1 tb- This suggests a 
breakdown of naive jet veto resummation at no later than N 4 LL. 

A corollary of the study is that a proper description of the Drell-Yan process with 
stringent jet vetoes must include entanglement between the two collinear sectors. It should 
be noted that our study only demonstrates the inevitability of entangling the two collinear 
sectors, while soft and ultra-soft gluons may still be allowed to factorize in some manner. 

Discrete symmetries play an important role in our approach. If factorization holds, 
the two collinear sectors are decoupled and have separate parity invariance, resulting in 
Z2 x Z-2 symmetry, which is richer than the Z2 global parity symmetry of QCD. Violation 
of factorization is revealed by the violation of extra discrete symmetries resulting from 
factorization. To construct a more realistic example that could be tested at colliders 
with unpolarized beams, one could exploit charge conjugation invariance: for example, the 
proton and the anti-proton have the same gluon beam function. If a future calculation 
shows that the Higgs production cross section (in the gluon fusion channel only) under 
a stringent jet veto in pp collisions is different, at leading power in C(py eto /Mhi ggs ), from 
the same quantity in pp collisions, it would be another manifestation of the violation of 
factorization. 
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A Definition of the beam function using QCD fields 

For completeness, we give a definition of the beam function using QCD fields. This ap¬ 
pendix is essentially a review and does not contain original work, because the definition 
agrees with the SCET definition [29] at least at low orders [35-37, 61]. The unsubtracted 
momentum-space beam function for a scalar parton (/> with Bjorken variable x\ and virtual¬ 
ity (ignoring transverse momentum components) oj \ = X\P + Mtr, for an incoming hadron 
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H\) with a large plus momentum component P + , 


^unsubtracted 




xP + 

2 


dwdw+ e -^ XlP+w + Mt rW +)/ 2( Hi ^ ( w +, w ~, 0 t ) 

/nr )h r ' ' 


V exp 


( 0 , 2 / , 0 T ) 


<K0)|Ri>, 


(A.l) 


where /r is the UV renormalization scale at which the matrix element is defined. This is 
directly analogous to Eq. (50) in [29] (with “+” and ” exchanged), except that the latter 
reference used the SCET collinear field with zero-bin subtraction. 

The eikonal beam function is defined by replacing the incoming hadron and the inter¬ 
polating field by a Wilson line in the “+” direction, 


5" (u l ,x 1 ,ri = 


xP + f dw dw + 


2n 2n 


—i(x\P+w +Mrfliu+)/2 


(0| w\ (w + ,w , Or) V exp 


dy . 


igA + (o, y , 0 T ) 


where we defined 


W\ (w + , w , wt ) = V exp 


dy + f , , _ \ 

/ -tt % 9 a \ w + y i w , w t) 

J —oo ^ 


Wi(0) |0>, 

(A.2) 

(A.3) 


Finally, we divide the Laplace transform of ^u nsubtracted with respect to w, B' 1 insubtracted , 
by the Laplace transform of L?® lkonal ) ^jeikonai^ to obtain the gauge-invariant moment-space 
beam function B\ (ui, x\, fi). The inverse Laplace transform of this result is the momentum- 
space beam function B\ (uj,x 
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